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Ion-dominated plasma and the origin of jets in quasars 
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Abstract. Low cooling plasmas associated with large kinetic energies are likely to be the origin of the kpc-extended 
and well collimated extra-galactic jets. 

It is proposed that jets are launched from a layer, governed by a highly diffusive, super-Keplerian rotating and 
thermally dominated by virial-hot and magnetized ion-plasma. The launching layer is located between the accre- 
tion disk and the corona surrounding the nucleus. The matter in the layer is causally connected to both the disk 
and to the central engine. Moreover we find that coronae, in the absence of heating from below, are dynamically 
unstable to thermal ion-conduction, and that accretion disks become intrinsically advection-dominated. 
We confirm the capability of this multi-layer model to form jets by carrying out 3D axisymmetric quasi-stationary 
MHD calculations with high spatial resolution, and taking into account turbulent and magnetic diffusion. The 
new multi-layer topology accommodates several previously proposed elements for jet-initiation, in particular the 
ion-torus, the magneto-centrifugal and the truncated disk - advective tori models. 
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1. Introduction 

Jets have been observed in many systems including active 
galaxies, X-ray binaries, black holes X-ray transients, su- 
persoft X-ray sources and young stellar objects (Konigl 
1997, Livio 1999, Mirabel 2001). Each of these systems is 
considered to contain an accretion disk, while jet-speeds 
have been verified to be of the order of the escape velocity 
at the vicinity of the central object (Mirabel 1999, Livio 
1999 and the references therein). Recent observations of 
the M87 galaxy reveal a significant jet-collimation already 
at 100 gravitational radii from the central engine, and that 
jet-launching should occur close to the last stable orbit 
(Biretta et al. 2002). 

Several scenarios have been suggested to uncover the 
mechanisms underlying jet-initiations and their connec- 
tion to accretion disks (Pudritz & Norman 1986). In most 
of these models magnetic fields (-MFs) are considered to 
play the major role in powering and collimating jets (e.g., 
the magneto-centrifugal acceleration model of Blandford 
& Payne 1982, the ion-torus model of Rees et al. 1982, X- 
point model of Shu et al. 1995, ADAF and ADIOS models 
of Narayan & Yi 1995 and Blandford & Begelman 1999). 

Previous radiative hydrodynamical studies without 
magnetic fields have confirmed the formation of a transi- 
tion layer (-TL) between the disk and the corona, governed 



Send offprint requests to: A. Hujeirat 



by thermally-induced outflows (Hujeirat & Camenzind 
2000). The aim of this paper is show that incorporating 
large scale magnetic fields (-MFs) manifests such forma- 
tion and dramatically strengthen the dynamic of the in- 
and out-flows. Moreover, the TL is shown to be an opti- 
mal runaway region where highly energetic ion-jets start 
off. The back reaction of jet-flows on the structure of the 
disk and on the corona surrounding the nucleus is inves- 
tigated also. 

The study is based on self-consistent 3D axi-symmetric 
quasi-stationary MHD calculations, taking into account 
magnetic and hydro-turbulent diffusion, and adopting the 
two-temperature description (Shapiro et al. 1976). This 
adaptation is fundamental as 1) the dynamical time scale 
around the last stable orbit may become shorter than the 
Coulomb-coupling time. Therefore turbulent dissipation, 
adiabatic or shock compression preferentially heat up the 
ions rather than electrons (7] oc p 2 J 3 , while T c cx p X J 3 )- 2) 
Taking into account that ions radiate inefficiently, having 
virial-heated ions in the vicinity of the last stable orbit is 
essential for the total energy-budget of large scale jets. 

2. Formation of the super-Keplerian layer 

Angular momentum in standard disks (Shakura & 
Sunyaev 1973, hereforth SS-disks) is transported out- 
wards mainly via small scale magneto-hydrodynamical 
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turbulence. Magnetic fields (-MFs) however were assumed 
to be weak, and remain so during the whole viscous evo- 
lution of the disk. Here we adopt a different view. Let 
ftot = ^hd + ^Mag = a(P gas + n B )/ri be a modified 
dynamical viscosity (see the caption of Figure 1 for elabo- 
ration). Unlike z^hd which transports angular momentum 
outwards, i-'Mag transports angular momentum in the 
vertical direction, provided that Bp is of large scale topol- 
ogy. Taking into account that i/hd in SS-disks decreases 
inwards, a transition radius i?Tr at which result^], the time 
scale of angular momentum removal by MFs at i?Tr is 
of the same order as that by the turbulent viscosity. Hence 



~a( — )( — ) ~1, 

Ttur r V A 



(1) 



where V s , Va arc the sound and transverse Alfven speeds 
(Va = \/ Bp Bt / in p) , respectively. Taking into account 
that H/r > a, we end up with Va > aV s - The last op- 
timistic inequality applies everywhere in the disk, irre- 
spective whether turbulence is mediated by 7Tb or not. 
Therefore, vertical transport of angular momentum via 
MFs is at least as efficient as a— viscosity. This is even 
more justified by the fact that Balbus-Hawley instability 
in Keplerian-disks amplifies initially weak fields to con- 
siderably large values, but remain still below equiparti- 
tion (Hawley et al. 1996). Imposing appropriate bound- 
ary conditions, and carrying the MHD-Box calculations 
with high spatial resolution, MFs could be amplified up to 
cquipartition (Ziegler 2002), yielding thereby T rcm < Ttur- 
Furthermore, the dynamo-action model proposed by Tout 
& Pringle 1992, if applied, would make r rcm even shorter. 

To be noted here that when taking a more realistic 
density and temperature stratification in global 3D MHD 
disk-calculations, vertical transport of angular momentum 
is inevitable (Arlt 2002). 

Efficient vertical transport of angular momentum rises 
the following important issues: 1) Turbulence in the disk 
need not be dissipated, or it might be even suppressed by 
the amplified Bp . This allows accretion to evolve without 
necessarily emitting the bulk of the their potential energy 
as radiation, and gives rise to energy re-distribution. 2) 
Accretion flows may not proceed as slowly as in SS-disk, 
but they may turn into advection-dominated^. This occurs 
because the time scale of angular momentum removal from 
the disk scales as: 



P V T H/BpB T 



r 3/2 



(2) 



where Vt(= rcos6* fl) is the angular velocity (Fig. 1). 
This indicates that angular momentum removal is more 
efficient at smaller radii. 



If t'Mag is too small, then magnetic fields are frozen-in to 
the gas approximately. In this case B ~ p 2//3 , and B must 
increase inwards, irrespectively, whether the disk is standard 
or advection-dominated 

2 The flow is said to be advection-dominated if transport via 
fluid-motion dominates viscous re-distribution. 



To maintain dynamical stability, angular momentum 
removal from the disk should be compensated by rapid 
advection from larger radii, i.e., r a( j v = T rem . Further, 
rapid and steady generation of in the disk yields 
Bp/B'Y = H&/r. The later two conditions imply that the 
radial velocity U T ~ Va , which means that the stronger the 
MF threading the disk, the more advection-dominated it 
becomes, and therefore the faster is the establishment of 
the super-Keplerian layer. 

Based on the present calculations (see the caption of 
Figure 1 and figures 2- 4), it is found that: 1) Angular ve- 
locity in the transition layer (TL) adopts approximately 
the profile Q ~ r~ 5 / 4 . 2) Energy dissipation is injected 
primarily into the ions that cool predominantly through 
fast outflows. 3) The generated toroidal magnetic field 
is quenched by a magnetic diffusion (reconnection) and 
fast outflows. The width of the TL is pre-dominantly 
determined through the transverse variation of the ion- 
pressure across the jet, i.e., f/\v = -Pi/VP; « 0.2 r, and so 
strongly dependent on whether the flow is a one- or two- 
temperature plasma. In the steady-state case, this implies: 

p~ r- 7/4 ,Ti - r- 1/2 ,U r ~ r~ 1/A ,B P /B T ~ const. (3) 

We note that U T adopts a profile and attains values 
similar to those in the innermost part of the disk. Provided 
that energy exchange between the matter in the disk and 
in the TL is efficient, the incoming matter can easily be 
re-directed into outwards-oriented motions. This implies 
that the Bernoulli number (Be) can change sign in dissi- 
pative flows. As Fig. 4 shows, Be is everywhere negative 
save the TL, where it attains large positive values, so that 
the ion-plasma can start its kpc-journey. Worth- noting is 
the resulting MF topology. Apparently, the outflow is suf- 
ficiently strong to shift the MF lines outwards, while the 
large diffusivity prevents the formation of large electric 
currents along the equator. In the corona however, MFs 
are too weak to halt the diffusive plasma in the dynam- 
ically unstable corona against gravity, and instead, they 
drift with the infalling gas inwards. In the case of very 
weak MFs (/3 < 0.1), our calculations indicate a consid- 
erably weak outflow. This is a consequence of the ten- 
dency of the MFs to establish a monopole like-topology, 
i.e, a one-dimensional MF-topology in which B# — > 0. In 
this case, the magnetic tension 7Tb becomes inefficient in 
feeding the matter in the TL with the angular momen- 
tum required for launching jets, indicating herewith that 
cold accretion disks alone are in-appropriate for initiating 
winds (Ogilvie & Livio 2001). 

Comparing the flux of matter in the wind to that 
in the disk, we find that Mw/ -Md = const, w 1/20. 
The angular momentum flux associated with the wind 
is Jw/jd = a,(M w / Aid) r 1 / 4 , where a is a constant of 
order unity. Consequently, at r = 3 x 10 2 i?sch, almost 
25% of the total accreted angular momentum in the disk 
re-appears in the wind. 

Why is the TL geometrically thin? 

In stratified dissipative flows the density scale height is 
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much smaller than the scale height of the angular velocity 
(Hujeirat & Camenzind 2000). Since the flow in TL 
rotates super-Keplerian, Coriolis forces act to compress 
the disk-matter and make its density scale height even 
smaller. This implies that advective-disks are geomet- 
rically thin, much thinner than what ADAF-solutions 
predict. 

On the other hand, unless there is a significant energy 
flux that heats up the plasma from below, as in the case 
of stars, heat conduction will always force the BH-coronae 
to collapse dynamically. To elaborate this point, let us 
compare the conduction time scale with the dynamical 
time scale along Bp-field at the last stable orbit of a 
SMBH: 

= Fro = 4-78 x 10 *p w T- io M, (4) 

r dyn KqT^ 

where pio, T;.io and A4$ are respectively in 10~ 10 gem -3 , 
10 10 K and in 10 8 A4® units. This is much less than unity 
for most reasonable values of density and temperature 
typical for AGN-environments. In writing Eq. 4 we have 
taken optimistically the upper limit c/y/3 for the veloc- 
ity, and set no = 3.2 x 10~ 8 for the ion-conduction coeffi- 
cient. When modifying the conduction operator to respect 
causality, we obtain T con d/i~dyn < U r /c, which is again 
smaller than unity. 

This agrees with our numerical calculations which rule out 
the possibility of outflows along the rotation axis, and in 
particular not from the highly unstable polar region of the 
BH, as ADAF-solutions predict. 

3. Summary 

We have presented a multi-layer model for initiating ion- 
jets in AGNs. The model agrees with a previous numerical 
study which confirmed the formation of thermally-induced 
outflows in the transition layer (Hujeirat & Camenzind 
2000). Here we have shown that incorporating the ef- 
fects of MFs manifests their formation and dramatically 
strengthen their dynamics. 

Three ingredients for initiating winds have been de- 
tected: 1) a highly diffusive plasma dominated by virial- 
hot ions, 2) large scale magnetic fields that efficiently 
transport angular momentum from the disk into the TL, 
where the plasma rotates super-Keplerian, and 3) an un- 
derlying advection-dominated accretion disk. 

Taking into account that the corona is dynamically 
unstable, adopting a large scale magnetic topology, and 
allowing ion-electron thermal decoupling appear to force 
accretion flows to undergo a global energy re-distribution: 
confined inflows (negative Bernoulli number) in the equa- 
torial region and in the corona, and thermally and 
magneto-centrifugally-driven outflows in the TL charac- 
terized through a positive Bernoulli number. This fea- 
ture may survive under different conditions: strong MFs 
suppress turbulence, weakening thereby the effect of the 
turbulent-viscosity and dominate the transport of angu- 
lar momentum. On the other hand weak MFs in rotat- 



ing stratified flows would be amplified via dynamo-action 
and reach equipartition, beyond which turbulence is again 
suppressed. This interplay between MFs and turbulent- 
viscosity, Balbus-Hawley and Parker instabilities may set- 
tle into an equilibrium state, in which inflows are si- 
multaneously associated with low-cooling out-flowing ion- 
plasma. 

We note that in the absence of thermal conduction and 
adopting the one-temperature description, low-viscosity 
radiativcly inefficient HD and MHD accretion flows be- 
come inevitably convection-dominated. Therefore, in the 
early phases of jet-initiation, CDAFs may play an impor- 
tant role in powering the jets in AGNs and microquasars 
(Abramowicz et al. 2002). 

The multi-layer model presented here accommodates 
some elements of BP82. In particular, we agree with BP82 
about the necessity for a super-Keplerian rotation of the 
plasma overlying the accretion disk. However, the plasma 
here is dominated by highly-diffusive and virial-hot 
ions; it does not require a special Bp— alignment with 
respect to the disk-normal to enable jet-launching, as 
ideal-MHD treatment requires. While our results agrees 
with the ion-torus model with respect to the necessity of 
2T-plasma to maintain the ions hot for a significant time 
of their propagation-life in the ISM, no signatures for the 
formation of ion-supported tori have been detected (Rees 
ct al. 1982). 

Our results differ from ADAF and ADIOs in several 
issues, and in particular with respect to 1) the existence 
of a layer adjusting to the disk, where the plasma is 
found to rotate super-Keplerian, 2) the configurations 
of the in- and the out-flows, 3) stability of the corona 
in the vicinity of the BH, 4) the transition from SS-disk 
to advection-dominated disks and 5) with respect to the 
essence of Bernoulli number in dissipative flows, i.e., a 
positive Bernoulli number is necessary but not sufficient 
for outflows (Abramowicz et al. 2000). 

Finally, we note that since the flow in the TL is highly 
dissipative (strengthen thermal and rotational coupling 
with the central nucleus), the innermost region of the disk 
rotates synchronously with Kerr black holes (due to the 
frame dragging effect), and since T rom decreases with ra- 
dius and depends inversely on Bt and Bp, we think that 
the plasma attached to the poloidal magnetic field would 
be forced to deposit its angular momentum to the plasma 
in the TL, thereby considerably enhancing the centrifugal 
power and ejecting the ion-plasma into space with rela- 
tivistic speeds. 
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Fig. 3. The distribution of the velocity field superposed on the logarithmic-scaled ratio of the ion- to electron- 
temperatures (red color corresponds to high ratios and blue to low-ratios). The lower figure is a zoom-in of the 
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Fig. 1. The model consists of 1O 8 M Schwarzschild BH 
at the center (its gravity is described in terms of the 
quasi-Newtonian potential of Paczynski & Wiita 1980), 
and an SS-disk (blue color, extending from r=l to r=20 
in units of the radius of the last stable orbit i.e., in 
3 x i?schwarzschiid, thickness — O.lx, an accretion rate 
of M = 0.01 x AlEdd, and a central disk temperature of 
T = 10~ 3 T V i r i a i at the outer radius). The ion-temperature 
Tj is set to be equal to the electron temperature T c 
initially.) The low-density hot corona (T = T V i r i a i, and 
density p(t = 0,r,9) = 10~ 4 p(i = 0,r,6» = 0) is set 
to envelope the disk. A large scale magnetic field is set 
to thread the disk and the overlying corona (blue lines, 



Vertical disribution 







1/4 at the outer radius, where P gas is 



the ion-pressure, P mag = B ■ B/8ir is the magnetic pres- 
sure, and B is the magnetic field whose components are 
(Bp, S T ) = (B1,B2,B T ).) The low-density hot corona 
(T = r virlal , and density p(t = O,r,0) = 10- i p(t = 
Q,r,8 = 0) is set to envelope the disk. A large scale 
magnetic field is set to thread the disk and the overly- 
ing corona (blue lines, (3 = P m ag/Pgas = 1/4 at the outer 
radius, where P gas is the ion-pressure, P mag = B ■ B/8n is 
the magnetic pressure, and B is the magnetic field whose 
components are (Bp, Bt) = (Bl, B2, Bt).) The numerical 
procedure is based on using the implicit solver IRMHD3 
to search steady-state solution for the 3D axi-symmetric 
two-tcmperaturc diffusive MHD equations in spherical ge- 
ometry (for further clarifications about the equations and 
the numerical method see Hujeirat & Rannacher 2001, 
Hujeirat & Camenzind 2000). The ion-pressure is used to 
describe the turbulent viscosity: ^turb = olP^/Q,, where 
a is the usual viscosity coefficient, and £1 is the angu- 
lar frequency. The magnetic diffusivity is taken to be 
equal to tWb- 250 x 80 strongly stretched finite volume 
cells in the radial and vertical direction have been used. 
Normal symmetry and anti-symmetry boundary condi- 
tions have been imposed along the equator and the ro- 
tation axis. Extrapolation has been adopted to fix down- 
stream values at the inner boundary. Non-dimensional 
formulation is adopted, using the reference scaling vari- 

3 ,f = 5 x 10 7 K, U = V S = 
Vsy/^np. The location of the 



ables: p = 2.5 x 10 12 gcm" 
jUgf/fr, - 1.23). B = 
transition layer (-TL), where the ion-dominated plasma is 
expected to rotate super-Keplerian and being accelerated 
into jets, is shown for clarity. 
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Fig. 2. The horizontal distribution of the normalized 
density p, angular velocity f2, radial and horizontal ve- 
locities Ur, V z , the radial and horizontal MF-components 
(51,52) = (B r , Be) = Bp and the toroidal MF- 
component Bt at r — 2.5. Note the density plateau, 
the positive radial velocity (outflows), the super-Keplerian 
rotation and the strongly enhanced strength of the MF- 
components in the TL. The inwards-oriented motions in 
the disk (inflows) and strongly increasing density towards 
the equator are obvious. 
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Fig. 4. A snap-shot of the distribution of the Bernoulli 
number in two-dimensional time-dependent calculations. 
The decrease from large to low positive values is repre- 
sented via yellow, green and red colors. The blue color cor- 
responds to negative values. The figure shows two ejected 
blobs of large positive energies in the TL. Ejection starts 
from the innermost region of the TL and evolves non- 
linearly. 



